Edge states versus diffusion in disordered graphene flakes 
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We study the localization properties of the wavefunctions at the Dirac point of graphene flakes 
with short range disorder which mixes the valleys. We accomplish this via the numerical calculation 
of the Inverse Participation Ratio(/P7?) and it's scaling which provides the fractal dimension D2. 
We show that the edge states which exist at the Dirac point of ballistic graphene (no disorder) 
with zig-zag edges survive in the presence of weak disorder with wavefunctions localized at the 
boundaries. We argue, that there is a strong interplay between the underlying destructive inter- 
ference mechanism of the honeycomb lattice of graphene leading to edge states and the diffusive 
interference mechanism introduced by the short-range disorder. This interplay results in a highly 
abnormal behavior, wavefunctions that are becoming progressively less localized as the disorder is 
increased, indicated by the decrease of the average {I PR) and the increase of D2. We verify, that 
this abnormal behavior is absent for graphene flakes with armchair edges which do not provide edge 
states. 



PACS numbers: 73.21. La, 73.22.Pr, 73.20.Fz, 73.22.-f 
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I. INTRODUCTION 

Grapheneii^ the first 2d metal ever made is an one 
atom-thick layer of carbon atoms arranged in a honey- 
comb lattice structure offering large practical advantages 
over conventional semiconductors, which make it an ex- 
cellent candidate for replacing silicon in future nanoelec- 
tronics. Most importantly graphene offers a vast field 
for fundamental theoretical work revealing phenomena 
such as, relativistic behavior of the electrons at the Fermi 
level(i?/ = 0), known as the Dirac point and topological 
phenomena like the edge states where the electron cur- 
rent flows from the sample boundaries like in topological 
insulators. 

Edge states have been studied in^"— through the theo- 
retical investigation of long stripes of graphene known as 
nanoribbons^^— while they have been experimentally ob- 
served in Ref. |Tl| through scanning tunneling microscopy 
(STM) and spectroscopy (STS) techniques. In Ref. |3, 
it is shown that nanoribbons with the so called zigzag 
edges exhibit zero energy edge states with wavefunctions 
concentrated at their borders. On the other hand, edge 
states are absent for the other possible type of edge mor- 
phology, armchair. So, zigzag nanoribbons exhibit vastly 
different electronic properties compared to the armchair 
nanoribbons, at the the Fermi level. In principle, the 
type of edge(zig-zag or armchair) at the boundaries of 
a graphene system is of crucial importance for its fun- 
damental electronic properties. This argument becomes 
also apparent when studying confined graphene struc- 
tures known as flakesi^rJ^ which have been experimen- 



tally studied In Ref. |]J|15| it is shown that 

graphene flakes exhibit edge states at the Fermi level de- 
pending on the edge type on their borders. For instance, 
trigonal and hexagonal flakes with zig-zag edges provide 
zero energy edge states which are absent for the corre- 
sponding shapes with armchair edges. 

Apart from the detailed edge morphology of graphene 
systems, another important factor that should be taken 



into account when studying their electronic properties is 
the presence of disorder which is an inevitable factor in 
any mesoscopic system. The main sources of disorder in 
graphene are the production method (synthesis) and the 
interaction with the supporting substrate j^ii^ Disorder 
can appear as lattice distortions like wrinkles, rippling, 
or impurities with various degrees of concentration com- 
ing from strains in the lattice or charge traps. Disorder 
is not always an undesirable factor , it can be useful also 
in applications, for example in spintronic devices through 
the interaction with the spi n^"'^^ which allows the manip- 
ulation of the magnetic properties of disordered graphene 
systems. Moreover, macroscopic graphene like structures 
(honeycomb) with controllable disorder have been shown 
to be achievable in Ref. i22l through microwave simula- 
tion of the electronic waves. The theoretical treatment 
of disordered graphene requires the introduction of dif- 
ferent models like short or long range disorder, simu- 
lating strong and weak concentration of impurities re- 
spectively. The type of disorder plays a crucial role on 
the localization properties of the wavefunctions. Long 
range disorder retains the separation between the two 
so called valleys, centered at the the two non equivalent 
Dirac points of pure graphene (E=0) at the corners of its 
hexagonal Brillouin zone, where the relativistic nature of 
electrons is revealed^ d^'^" . The separation between the 
two valleys in disordered graphene results in many in- 
teresting phenomena like anti-localization^ or minimum 
conductance^ir— . On the other hand, short range disor- 
der mixes the two valleys instead resulting in Anderson 
localization in Ref H^. In Ref. [33l - l36l a rough estima- 
tion of the localization properties of the wavefunctions 
can be derived through the study of energy level statis- 
tics of graphene flakes, verifying for example Anderson 
localization for short range disorder. So, there has been 
an extensive study of Anderson localization phenomena 
in graphene. 

However, the localization properties of the wavefunc- 
tions in the diffusive regime^^i*^ has been much less stud- 
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ied in graphene, especially concerning the cases where 
the edge states are present. The diffusive regime is de- 
fined in systems with short range disorder when the sys- 
tem's length is smaller than the localization length, where 
diffusive interference effects are known to dominate the 
behavior of the wavef unctions. This regime has been 
shown to hold phenomena that are absent in the local- 
ized regime (localization length larger than the system 
size), like multifractality of the wavefunctions^i. Since, 
these phenomena own their existence on the finite sys- 
tem size, it is reasonable to use confined structures like 
flakes for their investigation. Specifically, for graphene it 
is intersting to investigate these phenomena when edge 
states are present. So, our main goal in this paper is 
to investigate the edge states at the presence of disorder 
in the diffusive regime, through the study of disordered 
graphene flakes that provide edge states at the zero disor- 
der limit. Our analysis involves the numerical calculation 
of the Inverse Participation Ratio {I PR) and it's scaling 
behavior which gives the fractal dimension D-2r^^— char- 
acterizing roughly the volume of a wavefunction. Both 
measures combined give a rough picture about the form 
of the wave functions. Our calculations show evidence of 
the interplay between two mechanisms: the interference 
mechanism of the honeycomb lattice of graphene, lead- 
ing to localized wavefunctions at the borders namely the 
edge states, and the interference effects leading to dif- 
fusion of the wavefunction in normal disordered systems. 
The interaction of these two interference mechanisms has 
a large impact on the localization properties of the wave- 
functions. When edge states are present, we observe a 
decrease of the average {I PR) and an increase of D2 with 
increasing disorder, implying that the wavefunctions be- 
come progressively less localized. This is a highly ab- 
normal behavior compared to the normal 2d disordered 
systems, where the wavefunctions become naturally more 
localized with increasing disorder. We verify that the 
normal behavior is reproduced for disordered graphene 
flakes without edge states. 

The remainder of the paper is organized as follows. In 
section II, we introduce our numerical model based on 
the tight-binding framework along with the short-range 
disorder. In section III and IV, the numerical results 
for the various graphene flake shapes in the presence of 
short-range disorder are presented. We study two kind of 
shapes, triangular and hexagonal, with zigzag and arm- 
chair edges respectively. We discuss our results and con- 
clude in section IV. 



lowing Hamiltonian 

^ = Y1 + Mc'^+i + 4+lC^), (1) 

i i 

where cjj and c„ are the creation and annihilation op- 
erators for spinless fermions, t is the hopping element 
(t=1.0eV) and the random on-site potential, following 
the box distribution P{e) = l/w, in the range [-w/2,w/2] 
with w denoting the strength of the disorder. This type 
of disorder mixes the two valleys of graphene resulting 
in inter-valley scattering in our problem. We focus our 
study on graphene flakes with specific edge shape that 
are known to exhibit edge states at the limit of zero dis- 
order. These include shapes studied in Ref. [13 like the 
trigonal and the hexagonal, with zig-zag edges. Also, we 
extend our study on the same shape t ype s but with arm- 
chair edges, investigated also in Ref. [ij where the edge 
states are absent. 

Diagonalizing the Hamiltonians for each graphene flake 
we derive the wavefunction amplitudes i at each lat- 
tice site i for energy E, which allow us to calculate the 
inverse participation ratio (IPR) given by IPR{E) = 
\'^i{E)\'^. The participation ratio gives the informa- 
tion about the degree in which every site of the lattice 
is participating in the wavefunction. With its inverse we 
can get a rough estimation of the localization proper- 
ties. In general, a wavefunction with a dense distributed 
amplitude for example in the diffusive regime, will give 
values of IPR close to zero while a localized wavefunc- 
tion will give IPR values close to one. In addition, the 
scaling of IPR provides the fractal dimension D2 of a 
wavefunction through the relation IPR ^ L~^^ and is 
a characteristic measure of the volume it roughly occu- 
pies. For instance, when analyzing a two dimensional 
system, = 2 when the wavefunction amplitude is 
equally distributed on the whole lattice while D2 — Q 
for a localized wavefunction. In general, D2 belongs to 
a spectrum of fractal dimensions characterizing objects 
known as multifractals.Non integer values of D2 imply 
fractal wavefunctions^"— that are not either extended 
nor localized. 

The value of D2 gives a good estimation of the vol- 
ume a wavefunction is occupying although it does not 
uniquely identify it's form. For this reason, we illus- 
trate the wave- function probability amplitude |\E'i(i?)p 
for each case that we are studying. With these three 
measures the IPR, D2 and the wavefunction probability, 
we can get a fairly complete picture of the wavefunctions 
in disordered graphene flakes. 



II. MODEL 

For our analysis, we use a tight binding model for 
graphene with first nearest neighbor hopping and short- 
range disorder simulated by a random on-site potential 
on each lattice site. The model is described by the fol- 



III. TRIANGULAR GRAPHENE FLAKES 

In this section, we examine the wavefunction localiza- 
tion properties of disordered graphene flakes of trigonal 
shape with two different type of edges, armchair and 
zigzag (Fig. 1), via the analysis of IPR and its scal- 
ing which gives the fractal dimension Z?2. One example 
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FIG. 1 


The flake shap 


es we 


study, triangular and hexagonal. 



a) Equilateral triangular graphene flake with zigzag edges and 

b) arnichair edges, c) Hexagonal flake with zigzag and d) arm- 
chair edges. We choose the side length L of each shape, mea- 
sured in units of the lattice constant a, to characterize the 
flake size (linear length scale). The flakes with zigzag edges 
provide zero energy edge states at the zero disorder limit in 
contrast to the flakes with armchair edges. 



of a trigonal graphene flake with zigzag edges can be seen 
in Fig. la. We characterize the size of the flake by the 
base length of the triangular shape L in units of the lat- 
tice constant a — 2A2A, for example length is L—6a in 
Fig. la. In the case of zero disorder the zig-zag triaiigle 
exhibits zero energy edge states as discussed in Ref. [ij, 
with wavefunctions concentrated on the zigzag edges (see 
Fig. 3a). In Fig. 2, we plot the average value {I PR) 
versus the energy E near the Fermi level(-E = 0) for a 
specific flake of size L = 44 consisting of 2113 atoms and 
for different disorder strengths w— 0.5,1,1.5. The num- 
ber of disorder realizations is 5000. Orange color in the 
background represents the individual values of IPR for 
w=0.5 used to obtain the corresponding curve of (IPR). 
The < IPR > for w=0 can be also seen as large indi- 
vidual green colored dots, separated by large gaps. It is 
clear that for low energies, approximately until E ^ O.l 
the IPR is in average decreasing with increasing disorder 
strength. This is also evident in the inset of Fig. 2 where 
we show the corresponding values of IPR. The points 
concentrate progressively lower with increasing disorder 
resulting in lower values of {IPR) as we have seen in the 
main figure. Keeping in mind the two trivial limits of 
IPR, IPR ~ for an extended wavefunction and IPR = 
1 for a localized one, the overall behavior of IPR implies 
that the low energy wavefunctions of the zigzag triangle 
with disorder are becoming in average less localized as the 
disorder is increased. This behavior is highly abnormal 
and is absent in normal disordered materials, for instance 
the square lattice or the chain, where the wavefunctions 
become naturally more localized with increasing disor- 



der, with the IPR increasing in average. We can also 
observe an abrupt change of IPR starting at E ^ 0.1 
resulting in large fluctuations until E ^ 0.17. This is 
especially evident for w=0.5, in both the curve of (IPR) 
and the individual values of /Pi?(orange color). Below 
E ^ 0.1, we can distinguish a whole area of points with 
values of IPR in average much higher than the values 
corresponding to energies in the interval £' ~ 0.17 — 0.25. 

The energy area i? — 0.1 consists mainly of wave- 
functions that have their amplitude concentrated at the 
edges of the trigonal flake as shown in Fig. 3b, 3c where 
the wavefunction probability is plotted for E ^ 0.07 and 
different strengths of disorder. Comparing the edge state 
in Fig. 3a for w=0 with Fig. 3b we can see that weak 
disorder (w=0.5) localizes the wavefunction along a ran- 
dom area on the border of the flake acting in this way 
as a pertubation on the zero disorder limit studied in 
Ref. [Ij, where the edge state amplitude spreads peri- 
odically along the whole border. So, we can say that 
the edge states survive for weak disorder in the sense 
that the amplitude remains mostly concentrated at the 
border of the flake. The abnormal behavior that we dis- 
tinguished through the analysis of IPR in the energy 
area E ^ — 0.1 can be understood by looking at Fig. 
3c and 3d where we plot the wavefunction probability 
for stronger disorder w=1.5 and w=5. In Fig. 3c the 
amplitude although still mostly localized along border, 
has started extending across it while it penetrates more 
the flake, resulting in less localized wavefunctions and in 
lower IPR values as seen in Fig. 2. For strong disorder 
(w=5) in Fig. 3d the wavefunction becomes clearly local- 
ized inside the flake instead of the edges. For sufficiently 
strong disorder, the interference mechanism that leads 
to edge states in graphene is completely destroyed by 
the diffusive interference mechanism of the short-range 
disorder while even larger disorder leads to Anderson 
localization^^ with the wavefunctions becoming localized 
inside the flake. For energies above E ~ 0.17 in Fig. 2 
the localization properties of the wavefunctions change 
drastically spreading periodically along the whole bor- 
der flake indicated by IPR obtaining much lower values 
than for E ^ — 0.1. The transition from edge states 
to these extended states creates the large fluctuations in 
the energy interval [0.1 - 0.17]. 

We now continue our analysis by studying the scaling 
of IPR from which we can derive the fractal dimension 
D2 providing a rough estimation of the wavefunction's 
volume. In Fig. 4, we show D2 for different strengths of 
disorder with the values of the average (IPR) for each 
size being over 5000 disorder realizations and energies 
inside the window [0,0.1], where the edge states lie ap- 
proximately according to our previous analysis. Since 
IPR ~ we plot ln{{IPR)) versus ln{L) in or- 

der to get the exact value of D2. We can observe that 
the slope of each curve representing D2 increases with 
increasing disorder strength, implying that the volume 
occupied by the corresponding wavefunctions increases 
also. Moreover, ln{{IPR)) and consequently (IPR) for 
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FIG. 2: (color online). The average value of the Inverse Par- 
ticipation Ratio {{I PR)) versus the energy E for a trigonal 
graphene flake with zigzag edges of size L = 44(2113 sites), for 
disorder strengths w=0.5, 1, 1.5 and 5000 realizations, along 
with the w=0 case(green dots). The orange points in the 
background represent the individual values of IPR for w=0.5 
while the cases for the other disorder strengths are shown in 
the inset. {IPR) decreases with increasing disorder, implying 
less localized wavefunctions, visible also in the inset, where 
IPR concentrates on lower values with increasing disorder. 
There is an abrupt change of {IPR) at i? ~ 0.1 evident also 
in the background(orange points), coming from the transition 
from edge states concentrated at low energies to extended for 
higher energies. 



a) b) 




FIG. 3: (color online). The wavefunction probability ampli- 
tude for a flake of size L=44(2113 sites) for different dis- 
order strengths w=0.5,1.5,5 at energy E ~ 0.07(one disorder 
realization) along with a wavefunction for w=0 at i5 = 10^''. 
(a)For w=0 the amplitude is concentrated on the edges of the 
flake. (b)For w=0.5 the wavefunction becomes localized along 
the border. (c)For w=1.5, the amplitude penetrates slightly 
the flake, despite being mostly concentrated on the border, 
showing also abrupt fluctuations. (d)The wavefunction for 
large disorder strength w=5 concentrates in the flake's bulk. 
There is no sign of the edge states in this case, the wavefunc- 
tion becomes localized inside the flake instead. 



each individual size averaged over the energies and real- 
izations decreases, in agreement with the results obtained 
in Fig. 2 for the curves of {IPR) versus E. Additionally, 
D2 has non integer values implying multifractality inside 
the chosen energy window, evident from the abrupt fluc- 
tuations of the amplitude in Fig. 3. The values of D2 be- 
low one are reasonable considering the wavefunctions for 
w=0.5, 1.5 in Fig. 3b, 3c which are mostly concentrated 
along the border of the trigonal flake extending slightly 
inside. The overall behavior of D2 versus the disorder 
strength for the graphene triangle with zigzag edges can 
be seen in Fig. 14. So, we have found that although 
weak disorder preserves the edge states in trigonal flakes 
with zig-zag edges, its increase results in a highly abnor- 
mal behavior, the edge states become progressively less 
localized, extending inside the flakes. 

We contrast this result to the behavior observed for 
trigonal flakes with armchair edges shown in Fig. lb, for 
which edge states are absent for zero disorder according 
to Ref. [ij. {IPR} versus the energy E for flake size 
L=44. 5(2106 atoms)and for different disorder strengths 
can be seen in Fig. 5 along with their respective distri- 
butions of IPR in the background (orange points-w=0.5) 
and inside the inset(brown points-w=1.0, grey points- 
w=1.5). Also, the < IPR > for w=0 can be seen as large 
green dots. Clearly, as the disorder is increased, {IPR) 
increases inside the whole energy window [E = — 0.26] 




ln(L) 

FIG. 4: (color online). The scaling of IPR for a trigonal zig- 
zag flake(L=44) for different disorder strengths, averaged over 
energies in the interval [0, 0.1] and over 5000 realizations of the 
disorder. The slope of ln{{IPR)) versus ln(L) gives the fractal 
dimension D2 which allows the estimation of the wavefunction 
volume. It is clearly increased with increasing disorder while 
the points for IPR for small disorder strengths lie above the 
corresponding points for larger disorder in agreement with 
the behavior demonstrated in Fig. 2 considering the average 
behavior of IPR. 
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FIG. 5: (IPR) versus the energy for a trigonal graphene 
flake with armchair edges of size L=44.5(2106sites) for dis- 
order strengths w=0.5, 1, 1.5 and 5000 realizations, along 
with the w=0 case(green dots). The orange points in the 
background are the individual values of IPR for w=0.5, with 
the other cases shown in the inset. IPR increases in average 
with increasing disorder, implying more localized wavefunc- 
tions, a behavior observed in normal disordered systems like 
the square lattice. The gaps appearing at the distribution 
of IPR for w=0.5 (orange points) come from the respective 
gaps present for w=0. As the disorder is increased the gaps 
are disappearing, as shown in the inset. 



FIG. 6: The wavefunction probability amplitude for an arm- 
chair trigonal flake of sizeL=44. 5(2106 sites) for strength of 
disorder w=1.5 at iJ ~ 0.07. The amplitude spreads ran- 
domly all over the flake, this is a common picture of a diffusive 
wavefunction encountered in the diffusive regime of normal 2d 
systems. 



functions becoming more localized with increasing disor- 
der. 



IV. HEXAGONAL GRAPHENE FLAKES 



in which we perform our calculations, evident also in the 
insets with the points of IPR concentrating progressively 
in higher values. The gaps at specific energies for w=0.5 
are a consequence of the respective gaps appearing in 
the energy spectrum for zero disorder as seen from the 
< IPR > for w=0. For the zig-zag triangle in Fig. 2 
these gaps are overlapped by the spreading of the edge 
states to higher than zero energies where they lie for zero 
disorder. Apart from these gaps in Fig. 5, IPR behaves 
smoothly with the energy in contrast to the trigonal flake 
(Fig. 2) where we observed two regions with vastly dif- 
ferent values of IPR. In Fig. 6, we show a characteristic 
example of a wavefunction lying inside the energy win- 
dow of Fig. 5 a.t E ^ 0.07. The amplitude fluctuates 
wildly, randomly spreading on the whole flake, a com- 
mon picture of a diffusive wavefunction. 

The results for the fractal dimension D2 can be seen 
in Fig. 7 for different strengths of disorder. The value of 
IPR for each size is calculated as an average over differ- 
ent realizations of the disorder and over the energy win- 
dow [0, 0.2]. In contrast to the trigonal flakes with zigzag 
edges, D2 progressively decreases with increasing disor- 
der, as for normal 2d disordered systems. The overall 
behavior of D2 versus the disorder strength can be seen 
in Fig. 14. So, the trigonal flakes with armchair edges 
at the presence of short range disorder do not exhibit 
the abnormal behavior we encountered in the disordered 
trigonal flakes with zigzag edges. Instead, they behave as 
normal disordered metals with their corresponding wave- 



We now extend our analysis on hexagonal graphene 
flakes with zigzag edges which provide edge states at 
the zero disorder limit(see Fig. 9a) as in the case of 
the corresponding trigonal flakes^. Again, we are in- 
terested in the effect of the disorder on the edge states 
obtained through the scaling analysis of IPR. The over- 
all behavior of IPR versus the energy can be seen in 
Fig. 8 for a hexagonal flake of size L=18 consisting of 
1944 sites. The behavior is similar to that of the zig- 
zag triangle. For low energies below E ~ 0.1 where the 
edge states(see Fig. 9b, 9c) are concentrated, < IPR > 
and IPR are decreasing with increasing disorder meaning 
that the wavefunctions are becoming progressively less lo- 
calized as indicated by the wavefunction form in Fig. 9c. 
IPR changes drastically above E ^ 0.1 where extended 
states start to appear, obtaining much lower values de- 
spite the fact that the transition from the edge states to 
the extended states is smoother in this case than it is 
for the zigzag triangle because of the denser energy spec- 
trum for w=0(green points). Sufficiently stronger disor- 
der localizes the wavefunction inside the flake(see Fig. 
9d). As in the case of the zigzag trigonal flakes, D2 is in- 
creased as we increase the disorder (wavefunction volume 
increases), while its non-integer values imply multifrac- 
tality in agreement with the abrupt fluctuations of the 
amplitude in Fig. 9. It is clear that for weak disorder 
the edge states survive also in the case of the hexagonal 
flakes with zigzag edges and result in less localized states 
as the disorder is increased. So, the abnormal behavior 
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ln(L) 

FIG. 7: The scaling of IPR for a trigonal armchair flake for 
different disorder strengths, averaged over the energy interval 
[0, 0.2] and over 5000 realizations. The calculated slope D2 
characterizing the wavefunction volume is clearly decreased 
with increasing disorder while the points along the curves in- 
crease their values, in agreement with the average behavior of 
IPR observed in Fig. 4. 




FIG. 8: {IPR) vs. E for a hexagonal flake with zigzag 
edges of size L=18(1944 sites) for different disorder strengths 
and 5000 realizations, along with the case w=0(green points) 
and the individual values of IPR in the background (orange 
points-w=0.5) and the inset. {IPR) decreases with increasing 
disorder in agreement with the behavior of the IPR in the 
inset, implying in overall less localized wavefunctions. Around 
E ~ 0.15 there is a transition from edge states with high IPR 
values to extended states with low IPR values. The overall 
behavior is qualitatively similar to that of the trigonal flake 
with zigzag edge. 




FIG. 9: The wavefunction probability amplitude for a 
flake of size L=18(1944 sites) for different disorder strengths 
w=0.5,1.5,5 at i? ~ 0.07 along with the w=0 case at £ 
0.0005 . (a)The amplitude for w=0 clearly extends along the 
border of the flake. (b)For w=0.5, the amplitude remains 
concentrated on the border although localized in a specific 
area. (c)For w=1.5, the amplitude clearly starts penetrating 
the flake. (d)For strong disorder w=5, the wavefunction be- 
comes localized in the bulk. The overall behavior is similar 
to the trigonal flakes with zigzag edges. 




ln(L) 

FIG. 10: The scaling of IPR for a hexagonal zig-zag flake for 
increasing disorder strengths averaged over the energy inter- 
val [0,0.1] and 5000 realizations. The slope D2 is increased 
with increasing disorder while the points for IPRfor small dis- 
order lie above the corresponding points for larger disorder, 
in agreement with Fig. 7. 



we pointed out for the trigonal flakes with zigzag edges 
applies also for the hexagonal flakes with zigzag edges. 
We conclude that this effect is independent of the over- 
all flake shape and is related with the existence of edge 
states, governed by the detailed edge structure on the 
borders of the graphene flakes. 



To finalize our study, we consider the case of hexagonal 
flakes with armchair edges shown in Fig. Id. In this case, 
there are no edge states in the limit of zero disorder, in 
contrast to the flakes with zigzag edges. In Fig. 11, we 
show {IPR) for a hexagonal flake with armchair edges of 
size L=20.2 consisting of 2382 sites for disorder strengths 
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FIG. 11: {IPR) for a hexagonal flake with armchair edges 
of size L=20. 2(2382 sites) for different disorder strengths and 
5000 realizations along with IPR in the background and the 
inset including the w=0 case. Overall, IPR increases with 
increasing disorder like in the case of the trigonal armchair 
flake. This is the behavior observed in the diffusive regime of 
normal disordered systems. 




FIG. 12: The wavefunction probability for a hexagonal arm- 
chair flake consisting of 2382 sites(L=20.2) for strength of 
disorder w=1.5 at energy E ~ 0.07. The amplitude is spread 
on the whole lattice, fluctuating randomly,a characteristic ex- 
ample of a diffusive wavefunction. 



w=0,0.5,1.5 and 5000 realizations along with the corre- 
sponding IPR values in the background and the inset. 
Overall IPR increases as the disorder is increased like 
in the case of the trigonal shape with armchair edges. 
This is compatible with the behavior observed in normal 
metals, as we have already pointed out. Also in Fig. 13, 
we observe that D2 is clearly decreased. In Fig. 12, the 
corresponding wavefunction has the characteristic diffu- 
sive form with randomly fluctuating amplitude covering 
the whole lattice, compatible with the non integer val- 
ues of I?2, close to two. So, the hexagonal flakes with 
armchair edges at the presence of short range disorder, 
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FIG. 13: The scaling of IPR for a hexagonal flake with arm- 
chair edges for disorder strengths w=0.5, 1, 1.5 averaged in 
the energy interval [0, 0.2] and 5000 realizations. D2 behaves 
conventionally, decreasing with increasing disorder. 



do not exhibit the abnormal behavior we encountered in 
flakes with zigzag edges which provide edge states. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have presented numerical results for 
graphene flakes with short-range disorder that show a 
highly abnormal behavior for the localization properties 
of the wavefunctions when edge states are present. We 
observe a decrease of the average {IPR) and an increase 
of the fractal dimension D2 with increasing disorder, im- 
plying that the wavefunctions become roughly less lo- 
calized as the disorder is increased. We argue that the 
underlying mechanism that causes this behavior is the 
dynamic interplay between the destructive interference 
mechanism that produces edge states (localized wavefunc- 
tions at the borders) and the diffusive interference mech- 
anism, known to prevail in 2d mesoscopic systems with 
short range disorder when the system's length is smaller 
than the localization length. We have verified this behav- 
ior through the study of trigonal and hexagonal graphene 
flakes with zigzag edges where edge states are present. 
The abnormal behavior is absent for flakes with armchair 
edges which do not result in edge states. For sufficiently 
strong disorder all cases give wavefunctions that are be- 
coming progressively localized at the bulk of the flakes 
leading to Anderson localization. 

We would like to stress out that our conclusion comes 
in contrast to the widely accepted belief that there are no 
new phenomena concerning the localization properties of 
graphene at the presence of short range disorder. This 
argument is based on the fact that with short range dis- 
order the two non equivalent valleys of graphene at the 
Dirac point become mixed, so that, the relativistic na- 
ture of electrons is lost producing standard localization 
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FIG. 14: The fractal dimension D2 vs. the disorder strength 
w for trigonal and hexagonal flakes with zigzag and armchair 
edges. For the flakes with zigzag edges and w>0.25, D2 is in- 
creased with increasing disorder in contrast to the flakes with 
armchair edges for which D2 behaves conventionally, decreas- 
ing with the disorder. This means that when edge states are 
present the volume of the corresponding wavefunctions is in 
average increasing, they become progressively less localized 
with increasing disorder as we have shown through the anal- 
ysis of IPR. At the limit of zero disorder (w=0) D2 = 1 for 
both the zigzag triangle and hexagon since both flakes exhibit 
zero energy edge states extended along the flake's border. For 
w=0-0.25 localization of these edge states results in a steep 
decrease of D2, followed by the abnormal behavior we de- 
scribed. 



phenomena like Anderson localization. This is true for 
sufficiently strong disorder or when the length of the sys- 
tem is larger than the locahzation length, namely when 
the system lies on the localized regime. However, this 
view in our opinion ignores the interference effects which 
lead to diffusion of the wavefunctions. When edge states 
are present the consideration of this mechanism is very 
important and leads to unexpected behavior, as we have 
shown in the present paper. In future studies, we intend 
to extend our analysis in order to include effects like the 
magnetization of the edge a^°i^^ in graphene systems with 
disorder or to investigate the connection with topological 
insulators which has been shown to carry resemblance to 
graphene^^iii, due to the edge states mimicking the topo- 
logical property of the electron current flowing from the 
boundary surfaces in these materials. 



VI. ACKNOWLEDGEMENTS 

We would like to thank D. Katsanos, V. A. Gopar and 
S. N. Evangelou for useful discussions and careful reading 
of the manuscript. We also acknowledge the computer re- 
sources and assistance provided by the Institute for Bio- 
computation and Physics of Complex Systems (BIFI) of 
the University of Zaragoza. 



^ Geim, A. K. and Novoselov, K. S, Nature Materials 6,183- 
191 (2007). 

2 K.S. Novoselov et al.. Science 306, 666 (2004). 
^ K. Nakada, M. Fujita G. Dresselhaus, M. S. Dresselhaus, 
Phys. Rev. B 54, 17954 (1996). 

* Wakabayashi, K., M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

^ A. Akhmerov and CWJ Beenakker, Phys. Rev. B 77, 
085423 (2008) 

^ Castro Neto, A. H., F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim , Rev. Mod. Phys. 81, 109 
(2009). 

X. Li, X. Wang, L. Zhang, S. Lee, and H. Dai, Science 319, 
1229 (2008). 

* Liying Jiao, Li Zhang, Xinran Wang, Georgi Diankov, 
Hongjie Dai, Nature 458 (7240): 877-80 (2009). 

^ G. Brumfiel Nature, 458 (2009), p. 390. 
^° Michael S. Fuhrer, Nature Materials,9, 61 1-612 (2010). 
" C. Tao, L. Jiao, O. V. Yazyev, Y.-C. Chen, J. Feng, X. 

Zhang, R. B. Capaz, J. M. Tour, A. Zettl, S. G. Louie, H. 

Dai, and M. F. Crommie, Nat. Phys. 7, 616 (2011). 

Yamamoto, T.; Noguchi, T.; Watanabe, K. Phys. Rev. B 

2006, 74 (12), 121409. 
" W. L. Wang, S. Meng, and E. Kaxiras, Nano Lett. 8, 241 

(2008) 

H. P. Heiskanen and M. Manninen, J. Akola, New J. Phys. 
10, 103015 (2008). 



Akola J, Heiskanen H P and Manninen M., Phys. Rev. B 
77, 193410 (2008). 

L. A. Ponomarenko et al.. Science 320, 356 (2008). 

J. Guttinger, C. Stampfer, S. Hellmller, F. Molitor, T. Ihn, 

and K. Ensslin, Appl. Phys. Lett. 93, 212102 (2008). 

S. Schnez, F. Molitor, C. Stampfer, J. Guettinger, I. Sho- 

rubalko, T. Ihn, and K. Ensslin, Appl. Phys. Lett. 94, 

012107 (2009). 

Mucciolo, E. R., and C. H. Lewenkopf, J. Phys.: Condens. 
Matter 22, 273201 (2010). 
^° O. Voznyy, A. D. Giiclii, P. Potasz, and P. Hawrylak, Phys. 
Rev. B 83, 165417 (2011). 

Gleg V. Yazyev 2010 Rep. Prog. Phys. 73 05650 

S. Barkhofen, M. Bellec, U. Kuhl, F. Mortessagne, 

arXiv:1205.5488vl [cond-mat.mes-hall]. 

Ando, T., and T. Nakanishi, 1998, J. Phys. Soc. Jpn. 67, 

2857. 

'^'^ T. Ando, NPG Asia Mater, 1 (1) (2009), pp. 1721. 

K. Wakabayashi, Y. Takane, and M. Sigrist, Phys. Rev. 
Lett. 99, 036601 (2007). 

K. Wakabayashi, Y. Takane, M. Yamamoto, and M. 
Sigrist, Carbon 47, 124 (2009)[CAS]. 
P. W. Anderson, Phys. Rev. 109, 1492 (1958). 
P. W. Anderson, D. J. Thouless, E. Abrahams, and D. S. 
Fisher, Phys. Rev. B 22, 3519 (1980). 
S. J. Xiong and Y. Xiong, Phys. Rev. B 76, 214204 (2007). 
^° A. Rycerz, J. Tworzydlo, and C. W. J. Beenakker, Nat. 



9 



Phys. 3, 172 (2007) [CAS]. 

V.I. Falko and K.B. Efctov, Europhys. Lett. 32, 627 (1995); 
Phys. Rev. B 52, 17413 (1995). 

Dietmar Weinmann, Sigmund Kohler, Gert-Ludwig In- 
gold, and Peter Hanggi Ann. Phys. (Leipzig) (1999), Spec. 
Issue, SI-277-SI-280. 

I. Amanatidis and S. N. Evangelou, Phys. Rev. B 79, 
205420 (2009). 

Liang Huang, Ying-Cheng Lai, and Celso Grebogi Phys. 
Rev. E 81, 055203 (2010). 

M. Wimmer, A. R. Akhmerov, and F. Guinea, Phys. Rev. 
B 82, 045409 (2010). 

Adam Rycerz, Phys. Rev. B 85, 245424 (2012). 

H. Hentschel and I. Procaccia, Physica (Amsterdam) 8D, 



435 (1983). 

S. N. Evangelou, J. Phys. A 23, L317 (1990). 
Claudio dc C. Chamon, Christopher Mudry, and Xiao- 
Gang Wen Phys. Rev. Lett. 77, 4194,1996. 
Castillo H E, de C Chamon C, Pradkin E, Goldbart P M 
and Mudry C Phys. Rev. B 56, 10668,(1997). 
A. Rodriguez, L. J. Vasquez, K. Slevin, R. A. Romer Phys. 
Rev. B 84, 134209 (2011). 

Liwei Jiang, Yisong Zheng, Cuishan Yi, Haidong Li, and 

Tianquan Lu Phys. Rev. B 80, 155454 (2009). 

J. Li, I. Martin, M. Biittiker, and A. F. Morpurgo, Nature 

Phys. 7, 38 (2011). 

Biittiker, M., 2009, Science 325, 278. 



